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Filtering and Control of High Speed Motor Current in a 
Flywheel Energy Storage System 

Barbara H. Kenny and Walter Santiago 
National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

The NASA Glenn Research Center has been developing technology to enable the use of 
high speed flywheel energy storage units in future spacecraft for the last several years. An 
integral part of the flywheel unit is the three phase motor/generator that is used to accelerate 
and decelerate the flywheel. The motor/generator voltage is supplied from a pulse width 
modulated (PWM) inverter operating from a fixed DC voltage supply. The motor current is 
regulated through a closed loop current control that commands the necessary voltage from 
the inverter to achieve the desired current. The current regulation loop is the innermost 
control loop of the overall flywheel system and, as a result, must be fast and accurate over 
the entire operating speed range (20,000 RPM to 60,000 RPM) of the flywheel. The voltage 
applied to the motor is a high frequency PWM version of the DC bus voltage that results in 
the commanded fundamental value plus higher order harmonics. Most of the harmonic 
content is at the switching frequency and above. The higher order harmonics cause a rapid 
change in voltage to be applied to the motor that can result in large voltage stresses across 
the motor windings. In addition, the high frequency content in the motor causes sensor noise 
in the magnetic bearings that leads to disturbances for the bearing control. To alleviate 
these problems, a filter is used to present a more sinusoidal voltage to the motor/generator. 
However, the filter adds additional dynamics and phase lag to the motor system that can 
interfere with the performance of the current regulator. This paper will discuss the tuning 
methodology and results for the motor/generator current regulator and the impact of the 
filter on the control. Results at speeds up to 50,000 rpm are presented. 


Nomenclature 

iq S q-axis current (torque producing) in the rotor reference frame, amps 

ids d-axis current in the rotor reference frame, amps 

Vqds Commanded motor phase voltage vector in rotor reference frame, volts 

V qd s Commanded motor phase voltage vector in stator reference frame, volts 

iqds Motor phase current vector in the stationary reference frame, amps 

iqds Motor phase current vector in the rotor reference frame, amps 

A 

0 r Estimated rotor angle, degrees 

A 

oo r Estimated rotor speed, rad/sec 

T em Motor/generator torque, N-m 

A, a f Motor/generator back EMF constant, volt- sec 

L Inductance, henries 

C Capacitance, farads 

R Resistance, ohms 

Ki Integrator gain constant 

K p Proportional gain constant 

s LaPlace operator 

V dc DC bus voltage, volts 

f Fundamental frequency, Hz 
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I. Introduction 

F LYWHEEL energy storage units are an attractive alternative to traditional batteries in space applications that 
require large numbers of charge / discharge cycles, wide operating temperature ranges, or peaking power for 
short periods. In addition, multiple flywheels can be used on a satellite to provide both energy storage and attitude 
control, thus combining two subsystems into one. Flywheels are composed of multiple subsystems including the 
high inertia flywheel rotor, the motor/generator, the magnetic bearings and the auxiliary bearings as shown in Fig. 1. 
The NASA Glenn Research Center has been developing advanced technologies for these subsystems over the last 
several years and has demonstrated energy storage at 60,000 rpm with one unit and combined single axis attitude 
control and energy storage using two units [1,2]. 



Housing 


Composite Rotor 


A uxiliary Bearings 


Magnetic Bearings 


Motor/Generator 


Figure 1. Flywheel Unit and Components 


One important area of research is the 
development of the motor/generator 
controls. Algorithms have been developed 
to control the motor/generator such that 
the flywheel can store energy in charge 
mode and supply energy to loads in 
discharge mode while regulating the DC 
bus voltage [1]. Additional algorithms 
have been developed to combine the 
attitude control and energy storage 
functions of two separate flywheel units 
[2,3]. The innermost loop of the flywheel 
control for all of the algorithms is a 
current loop that is based on the field 
orientation control technique in which 
current is directly proportional to motor 
torque [4,5]. A block diagram of the field 
orientation control of the flywheel 
motor/generator is shown in Fig. 2. 

In the field orientation technique, a 


mathematical transformation is performed on the current variables such that all of the control is done in a reference 
frame that is synchronous with the rotor. This allows the control quantities to become DC values and thus the 
bandwidth of the current regulator is not affected by the fundamental frequency of the machine. In addition, if the 
rotor reference frame d-axis current, id s , is regulated to zero, the machine torque becomes linearly proportional to the 
rotor reference frame q-axis current, iq S as shown in Eq. (1) where L af is the back EMF constant of the machine. Thus 
the machine torque response is directly proportional to the current response. 


"Cm 



( 1 ) 


To perform the reference frame transformation properly, continuous rotor angular position feedback information 


is necessary. In typical industrial implementations, this information is provided by an additional sensor such as an 


encoder or a resolver. For the flywheel 
application, due to the high operating 
speed required, there is no position or 
speed sensor on the shaft. A “sensorless” 
position and speed estimation algorithm 
provides the necessary information by 
using the measured motor phase current 
values and the commanded motor voltage, 
V q ds, values to continuously estimate rotor 
position and speed. This is described in 
detail in reference [6]. 

As the innermost loop of the entire 
flywheel control system, the algorithm 
shown in Fig. 2 must deliver a fast and 
accurate torque response over the entire 



Figure 2. Motor Control Block Diagram 
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operating speed range of the motor/generator. The key component to accomplishing this is the current regulator 
because the current is directly proportional to the torque as shown in Eq. (1). This paper discusses the current 
regulator performance with respect to the effect of the AC filter, proper tuning, and the operating speed with results 
presented for 20 kRPM and 50 kRPM operation (60,000 RPM was not demonstrated for this paper but has been 
documented previously [7,8]). 

II. AC Filter 


The AC filter is necessary in this application for a number of reasons. First of all, it reduces the dv/dt stress that 
the motor winding sees due to the high frequency (65 kHz) switching. Switching at a high frequency is necessary to 

create a high fidelity 1 kHz fundamental frequency at 
the full speed value of 60,000 RPM. Secondly, the 
filter reduces the current ripple in the motor. The 
current ripple is significant without the filter because 
the machine inductance is relatively low. With the 
filter, this ripple is approximately 1 amp, peak to 
Figure 3. Filter Circuit Diagram peak. It also reduces the common mode noise in both 

the motor and the magnetic bearings due to the 
parasitic capacitive coupling between the windings 
and the motor housing. The eddy current sensors 
used to sense the rotor shaft position for the magnetic 
bearings are particularly sensitive to high frequency 
noise. Further detailed discussion of the design and 
benefits of the AC filter can be found in [9]. 

The AC filter itself consists of two cascaded 
filters: a 65 kHz “trap” filter tuned to reduce the 
harmonics at the switching frequency, and a four 
Table 1. Filter Component Values pole, two stage low pass filter tuned with a 10 kHz 

corner frequency. The circuit diagram of one phase 
of the three phase filter is shown in Fig. 3. V A is the voltage at the inverter terminal, V a is the voltage at the motor 
terminal and n is the filter neutral. The motor neutral is not connected but the filter neutral is connected to the 
midpoint of the DC bus capacitor on the DC side of the inverter. The filter has the component values given in Table 
1. 




III. Synchronous Frame Current Regulator 

The synchronous frame current regulator [10] shown in the block diagram in Fig. 2 consists of two PI 
controllers, each one acting on an axis of current, as shown in Figs. 4 and 5. The outputs of the current regulators 
result in a rotor reference frame voltage command that is converted to a phase voltage (line-neutral) command in the 


stationary reference frame. The conversion is accomplished using the estimated rotor angle, 0 r , as shown in Fig. 2. 


The commanded phase voltage is synthesized from the DC voltage by the high 
frequency pulse width modulation (PWM) switching of the inverter. 

From a controls perspective, the simplest model of the motor is an R-F 
circuit in the rotor reference frame as described by Eqs. (2) and (3). 


di i = J_( 
dt F q V 
dids = J_ ( 
dt F d V 


Vq S - R s iq S 
Vds - R s ids 


( 2 ) 

( 3 ) 


Using FaPlace variables and rearranging terms results in expressions for the 
transfer functions as shown in Eqs. (4) and (5) . 

1 



( 4 ) 



Figure 4. q-Axis Synchronous 
Frame Current Regulator 



Figure 5. d-Axis Synchronous 
Frame Current Regulator 
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motor model 
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( 5 ) 


Figure 6. q-Axis Controller and Motor Block Diagram 


controller model 


motor model 


The motor model can now be combined 
with the current regulator model into one 
block diagram so that the initial tuning of the 
regulator can be found as shown in Figs. (6) 
and (7). The PI current regulator has been 
rewritten in a form to show the zero of the 
controller and the inverter is modeled with a 
transfer function of 1. Note that with the 
simple R-L model of the motor in the rotor 
reference frame, there is no coupling between 
the d- and q- axes. 

The simple inverter model with a transfer 
function equal to 1 is reasonably accurate for the fundamental frequencies as long as the magnitude of the 
commanded voltage does not exceed a maximum value that is dependent on the DC bus voltage magnitude and the 
modulation technique. The magnitude of the commanded voltage is given in Eq. (6) and is the same in both the 
rotating and stationary reference frames. The modulation technique used in the flywheel system is known as space 
vector modulation. For this modulation technique, the maximum commanded voltage magnitude that allows the 
inverter to remain in the linear region of control (transfer function equal to 1) is given in Eq. (7) [11]. 



Figure 7. d-Axis Controller and Motor Block Diagram 


|V£ s |=V(Vq;) 2 + (V£) 2 


|v&| =V (<) 2 + (vf s ) 2 


|Vqds| — Vmax 


Vdc 

V 3 


( 6 ) 

(V) 


Based on Figs. 6 and 7, the transfer function of the q- and d-axis currents are given in (8) and (9), respectively. 


li 

ids 



( 8 ) 


( 9 ) 


For the purpose of tuning the current regulator, L s , the average of the d- and q-axis inductances, is used. With this in 
mind, if the gains in Eqs. (8) and (9) are set according to Eq. (10), then the transfer functions reduce to the values 
shown in Eqs. (11) and (12). 


Rs 

L s 


Kl 

K P 


ids 

i£ 



( 10 ) 

( 11 ) 


( 12 ) 


The bandwidth of the response can then be set according to Eq. (13) with the selection of K p which sets the poles 
of Eqs. (11) and (12) to the desired location. The bandwidth in this equation corresponds to the response of the 
motor torque and is theoretically independent of the fundamental frequency and operating speed of the machine. 
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f bw =T^ Hz (13) 

Generally, the phase resistance and inductance of the motor are used in Eqs. (10) through (13) to calculate the 
current regulator gains and bandwidth. Even if perfect pole-zero cancellation is not achieved in Eqs. (8) and (9), 
either the response is adequate for a particular application, or the gains can be tuned heuristically starting from the 
calculated values until the desired response is achieved. However, if the filter introduces a large impedance into the 
circuit, the response can be significantly different than expected. The initial tuning can be improved by modifying 
the inductance and resistance values used in Eqs. (10) through (13) to include the effect of the filter impedance. 

The filter impedance can be calculated using a Thevenin equivalent impedance approach and the circuit diagram 
of Fig. 3. The Thevenin equivalent voltage is basically the fundamental voltage from the inverter with significantly 
reduced 65 kHz switching content and harmonics. It has very little phase lag or magnitude change within the 
operating frequency range (0-1000 Hz) as discussed in [9]. The Thevenin equivalent impedance of the filter shown 
in Fig. 3 is a 6 th order equation that is given in the Appendix. At the frequencies of interest, however, it is 
predominantly inductance, and this is most easily seen by actually measuring the impedance of the circuit using an 
impedance analyzer and then calculating the equivalent inductance and resistance. Three possible circuit 
configurations were measured: the motor only, the motor with the two stage filter of Fig. 3 but without the L-C trap 
portion, and the motor with the entire filter of Fig. 3. The results are given in Table 2. 


Frequency 

(Hz) 

Motor Only 

Motor plus filter without 
trap 

Motor plus filter with trap 
(Fig. 3) 

Req (mQ) 

L e q (pH) 

Req (mQ) 

L e q (pH) 

Req (mQ) 

L e q (pH) 

100 

33 

39 

66 

64 

84 

140 

200 

36 

38 

68 

63 

87 

139 

300 

38 

37 

71 

62 

90 

138 

400 

41 

37 

74 

62 

94 

138 

500 

44 

36 

78 

62 

99 

137 

600 

47 

36 

82 

62 

104 

137 

700 

50 

35 

87 

62 

109 

137 

800 

54 

35 

93 

62 

116 

137 

900 

57 

34 

99 

62 

122 

137 

1000 

60 

34 

106 

62 

130 

137 

Average 

46 

36 

82 

62 

104 

138 


Table 2. Circuit Resistance and Inductance Values Calculated from Impedance Measurements. 

Based on the measured resistance and inductance, two sets of gains can be calculated for various bandwidths 
using (10) and (13) as shown in Table 3. 


Bandwidth 

Motor Only Parameters 

Motor plus two stage filter (no 
L-C trap) 

Motor plus two stage & 
trap filter (Fig. 3) 


Kp 

Ki 

Kp 

Ki 

Kp 

Ki 

1 kHz 

.2 

544 

.39 

515 

.87 

654 

1.5 kHz 

.3 

820 

.58 

773 

1.3 

976 

2 kHz 

.4 

1088 

.78 

1030 

1.73 

1308 


Table 3. Current Regulator Gains for Various Bandwidths. 


IV. Experimental Results 

For all of the experimental results, the two stage with F-C trap filter (Fig. 3) was used between the output of the 
inverter and the motor as shown in Fig. 2. However, the first two sets of results presented herein use current 
regulator gains based on the motor parameters alone, without consideration of the effect of the filter. The 
performance consequences of not including the filter impedance in calculating the gains can then be seen. The rest 
of the results use the gains based on the entire circuit impedance measurement: motor plus two stage and F-C trap 
filter (the last two columns of Table 3). All of the results use the gains calculated for a 2 kHz bandwidth. 

Two speeds within the operating range of the flywheel system, 20,000 RPM and 50,000 RPM, were tested. In 
general, either set of gains worked well at 20,000 RPM. The response was fast and there was very little overshoot or 
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time, seconds 


Figure 8: Motor Phase Current Response to Step 
Change in q-Axis Current Command at 20 kRPM. 


oscillation. In contrast, operation at 50,000 RPM was 
more challenging and required modifications of the 
controller to achieve the desired performance. To best 
demonstrate the performance improvements resulting 
from the control modifications, the subsequent data 
presented is at 50,000 RPM. The control modifications 
also improve performance at low speeds but the 
improvement is more evident at higher speeds and thus 
the high speed data is presented herein. 

Figures 8 through 12* show 20,000 RPM operation 
using the motor alone gain parameters. Figure 8 shows 
the phase current response to a step change command in 
motor q-axis current from 1.5 amps to 20 amps. 
Figures 9 and 11 show the motor q-axis current (and 
estimated torque) and the motor speed. The motor 
accelerates at a constant rate after the step change, 
indicating a constant motor torque as shown on Fig. 11. 
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ns 
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0.0 
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Figure 9. Motor Commanded (red) and Actual (blue) 
q-axis Current and Equivalent Torque at 20 kRPM. 



Figure 11. Motor Speed for Step Change in 
Commanded Current at Time t=0 at 20 kRPM. 



Figure 10. Expanded Scale of Motor q-axis Current 
Response at 20 kRPM. 



Figure 12. Expanded Scale of Motor d-axis Current 
Response at 20 kRPM. 


In all of the plots in this paper, the red trace is the commanded value and the blue trace is the measured value. 
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time, seconds 

Figure 13. Motor Phase Current Response to Step 
Change in q-Axis Current Command at 50 kRPM. 


Figures 10 and 12 show an expanded time scale 
version of the q- and d-axis current response. It can be 
seen that the q-axis current rises to the commanded 
level with a response time on the order of 2 kHz but 
with some overshoot. The d-axis current exhibits a 
disturbance when the q-axis step change occurs. This 
is due to coupling between the d- and q-axes that is not 
accounted for in the simple R-L model of the motor 
and will be discussed in Section V. 

Figures 13 through 17 give the results at 50,000 
RPM, again using the gains based only on the motor 
parameters. This corresponds to an 833 Hz 
fundamental frequency of current and voltage. It can 
be seen that at this speed there is significant oscillation 
as the current reaches the commanded value. Figure 
15 shows that the overshoot of the q-axis current 
reaches almost 150% of the commanded value and the 
oscillation takes over 20 msec, to dampen out. The 



Figure 14. Motor Commanded and Actual q-axis 
Current and Equivalent Torque at 50 kRPM. 



Commanded Current at Time t=0 at 50 kRPM. 


Figure 15. Expanded Scale of Motor q-axis Current 
Response at 50 kRPM. 



Figure 17. Expanded Scale of Motor d-axis Current 
Response at 50 kRPM. 


NAS A/TM— 2004-2 13343 


7 




time, seconds 

Figure 18. Motor Phase Current Response at 50 
kRPM, with Filter Gains. 


current regulator performance using the motor 
parameters only gains is not acceptable at high speeds 
and may even become unstable for larger disturbances. 
In contrast, Figs. 18 through 22 show the current 
regulator performance at 50,000 RPM using the gains 
calculated with the total measured impedance (motor 
plus the two stage filter with L-C trap, Fig. 3, and a 2 
kHz bandwidth from Table 3). The overshoot is still 
present but it is reduced in magnitude and there is no 
period of extended oscillation. The q-axis current 
settles to the commanded value in less than 5 msec, 
and the rise time is very fast, achieving the desired 2 
kHz bandwidth. However, the d-axis current still 
exhibits a large transient during the step change in the 
q-axis command. This transient can be reduced by 
augmenting the current regulator to compensate for the 
cross coupling between the q- and d-axes that is not 
represented in the simple models of Eqs. (4) and (5) 


35 
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Figure 19. Motor Commanded and Actual q-axis 
Current and Equivalent Torque at 50 kRPM with 
Filter Gains. 



Figure 21. Motor Speed during transient at 50 kRPM 
with Filter Gains. 



time, seconds 

Figure 20. Expanded Motor q-axis Current Response 
at 50 kRPM with Filter Gains. 



Figure 22. Expanded Motor d-axis Current Response 
at 50 kRPM with Filter Gains. 
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Figs. 6 and 7. The cross-coupling causes the d-axis current to exhibit a response to a q-axis change and the q-axis 
current to exhibit a response to a d-axis change. For a perfectly decoupled system, a change in the q-axis current 
should cause no change in the d-axis. 


V. Back EMF Decoupling 


The current response can be improved by designing the current regulator using the complete motor model 
instead of the simple R-L model previously given in Eqs. (4) and (5). The complete state equations are given in (14) 
and (15) for the permanent magnet machine modeled in the rotor reference frame [5]. 


dt 

dids 

dt 



Vq S - Rsiqs “ LdO) r ids “ CX) r A, 


Vds “ Rsids Lq(O r iq, 


(14) 

(15) 


In addition to the resistance and inductance terms of the simple R-L model, there is now a cross-coupling term in 
each equation that is dependent on the inductance, the speed and the current in the other axis. In addition, there is 
the back EMF term in the q-axis equation, oo r X a f, that represents the voltage produced by the rotor magnets as they 
rotate through the stator windings. 

The effect of the cross-coupling can be removed by adding terms to the current regulator that anticipate and 
approximately cancel the coupling terms using a technique known as back EMF decoupling [5]. If the motor 
parameters are known exactly, the result of back EMF decoupling is to reduce the motor state equations to the 
simple R-L model given earlier in Eqs. (4) and (5). This is shown in (16) and (17) where the A symbol indicates the 
estimated value of the parameter. The block diagram of the total system, current regulator and motor model, is 
shown in Fig. 23. The right half of the block diagram is the motor model and the left half is the current regulator 
with back EMF decoupling. The gain values, K p and K i? remain the same as given in Table 2. 


di, 

dt 


il 1 ( r . r . r ) 1 ( A . r A \ if r 

T \ Vq S - Rsiqs " LdO) r ids - Q) r A. a f / A V Ld(O r ids + (X) T 'h a f/ ~ T \ Vq S - R s iqs 


dids = J_( 

dt L d V 


Vds - Rsids + LqCOj 


r \ 1 ( A \ 1 ( 

'dqs / “ a ' Lq(O r iqs / ~ y V 

r -*^d 


q UJ r 1 qs 


V d r s - R s id r s 


(16) 

(17) 



Figure 23. Back EMF Current Regulator and PM Motor Modeled in the Rotor Reference Frame. 
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time, seconds 


Figure 24. Motor Phase Current Response at 50 
kRPM, with Filter Gains and Back EMF Decoupling. 


The results using the controller with the back EMF 
decoupling added are shown in Figs. 24 through 28. 
The flywheel speed is 50,000 rpm and the gains are the 
same as for the results shown in Figs. 18 through 22. 
The decoupling current regulator reduces the d-axis 
overshoot as can be seen by comparing Fig. 28 to Fig. 
22. In addition, the q-axis current peak is also reduced 
as shown in Fig. 26 versus Fig. 20. However, Fig. 24 
shows a greater phase current overshoot than in Fig. 18 
in spite of the improvements seen in the d- and q-axes. 
In the ideal case, an improvement in the q- and d-axis 
response should lead to a corresponding improvement 
in the phase current response. 

The most likely explanation for this discrepancy is 
that there is an error between the actual and estimated 
rotor angle, 0 r , that is used in the reference frame 
transformations necessary for the field orientation 


35 


30- 

25- 


0.45 



time, seconds 


Figure 25. Motor Commanded and Actual q-axis 
Current and Equivalent Torque at 50 kRPM with 
Filter Gains and Back EMF Decoupling. 



Figure 27. Motor Speed during transient at 50 kRPM 
with Filter Gains and Back EMF Decoupling. 



time, seconds 


Figure 26. Expanded Motor q-axis Current Response 
at 50 kRPM with Filter Gains and Back EMF 
Decoupling. 
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Figure 28. Expanded Motor d-axis Current Response 
at 50 kRPM with Filter Gains and Back EMF 
Decoupling. 
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control technique. The estimated rotor angle is based primarily on the assumption that the commanded voltage, V q d s , 
is equal to the actual motor phase voltage [6]. However, the filter adds additional impedance between the inverter 
and the motor as already described. This impedance has an associated phase lag that adds to the controller sampling 
and processing time lag. The total phase shift will appear directly as an error in the estimation of the rotor position 
angle, 0 r . The result of using an inaccurate rotor angle is to introduce additional coupling between the d- and q-axis 
in the torque response of the machine and also to deviate from the most efficient (measured as maximum torque per 
amp) operating point. Identifying and quantifying the impact of rotor angle estimation errors is an area of future 
work. 


VI. DC Bus Voltage Limit 

As discussed in Section III, the motor control is based on the assumption that the transfer function between the 
commanded voltage and the output fundamental voltage from the inverter is equal to 1 . This is true as long as the 
commanded voltage is less than the maximum value given in Eq. (7). In the flywheel system with a DC bus voltage 
of 125 volts, the maximum peak value is 72 volts, line to neutral. This peak value ultimately limits the speed of the 
machine because as the speed increases, so does the required motor phase voltage. This can be seen from the motor 
state equations given in Eqs. (14) and (15) rearranged to see the necessary voltage for various operating conditions 
as shown in Eqs. (18) and (19). Equation (20) describes the peak phase voltage magnitude. 



Figure 29. Commanded Phase Voltage During Step 
Change in Current Command at 20 kRPM. 
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dJi 

dt 


(18) 

(19) 

( 20 ) 


For steady state operation (constant speed), the 
derivative terms in (18) and (19) become zero. The 
commanded d-axis current is zero and if the rotor 
angle used in the reference frame transformation is 
correct, the actual d-axis current will also be zero and 
thus the i^ terms will be equal to zero in the equations. 
The voltage equations then reduce to Eqs. (21) and 
( 22 ). 

( 21 ) 


V q s F-siqs ^r^af 

V ds L q O) r iqs 


( 22 ) 


It can be seen from Eqs. (21) and (22) that at 
steady state, the magnitude of the voltage applied to 
the machine will be dominated by the back EMF 
voltage (the oo r X a f term). This is because the voltage 
drop across the motor resistance in Eq. (21) and the 
voltage drop across the inductance in Eq. (22) are 
small compared to the back EMF voltage, particularly 
at high speeds. During a step change in the current 
command, however, the voltage magnitude will 
increase sharply due to the non-zero value of the 
derivative term in Eq. (18). After the step, during the 
acceleration of the machine, the voltage magnitude 
will continue to be higher than in steady state because 
the current, i qs , increases to provide the torque 
necessary for acceleration. 


Figure 30. Commanded Phase Voltage During Step 
Change in Current Command at 50 kRPM. 

The voltage magnitude during the steady state, the step change transient, and the acceleration phases can be seen 
in Figs. 29 and 30 for 20,000 RPM and 50,000 RPM operation, respectively. Both plots are for a 2 kHz bandwidth 
system using back EMF decoupling and the gains given in Table 3 (with the filter of Fig. 3). In addition, Fig. 30 
corresponds to the same test run as reported in Figs. 24 to 28. It can be seen from the figures that the voltage spike 
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due to the step change in the current command at time t=0 is within the maximum available voltage for 20,000 RPM 
operation but exceeds it for approximately a millisecond for 50,000 RPM operation. During this time, the inverter 
provides the maximum voltage possible but cannot provide the commanded value. However, the current is able to 
eventually reach the correct operating point even though the commanded voltage is not available for the first 
millisecond of the transient. 

It can also be seen from the figures that for the 20,000 RPM operation case, the voltage magnitude behaves as 
described above: close to the back EMF voltage at steady state (t < 0) and slightly higher than the back EMF voltage 
during acceleration (t > 0). For the 50,000 RPM case, the voltage magnitude is approximately equal to the back 
EMF voltage in steady state but it clearly decreases during acceleration. This would seem to contradict the 
explanation given above, however, that derivation is based on the assumption that the reference frame 
transformations are accomplished correctly through accurate knowledge of the rotor angle. In the present 
implementation of the controller, a phase lag between the commanded voltage and the voltage applied to the motor 
would cause a rotor angle estimate that lags the actual angle and results in a non-zero (negative) id s term. A negative 
ids will reduce the voltage magnitude as can be seen from Eqs. (18), (19) and (20). This is the most likely cause of 
the voltage magnitude reduction observed in Fig. 30. As stated previously, work is presently underway to improve 
the rotor angle estimate and quantify any errors, particularly at high speeds, and it is expected that a more accurate 
estimate will resolve the discrepancies between the expected and measured behavior. It is interesting to note, 
however, that the controller still continues to provide an excellent overall response in spite of the suspected rotor 
angle estimation error. 

Returning to the problem of the commanded voltage exceeding the available voltage, one possible solution is to 
limit the value of the derivative term in Eq. (18) by limiting the rate of change of the commanded current. If a 2 
kHz response is desired, it is not necessary to change the current command in a step fashion but rather limit it to a 
ramp function that represents a command change within 2 kHz. Figures 31 through 36 show the results for a 2 kHz 
bandwidth response, tuned with the gains of Table 3, and a 60 kA/sec ramp for the commanded value of current, iq S , 


as it transitions from 1.5 to 20 amps. 


This means that the derivative term in Eq. (18), 



, will be approximately 


equal to 60 kA/sec. Figure 31 shows the peak voltage command to be significantly less than the case without the 
ramp command and now well below the maximum available voltage. In addition, the current response is improved. 
There is less overshoot on the d-axis current (Fig. 36) and the q-axis current shown in Fig. 34 still has a 2 kHz 
response with an approximate settling time of 2 msec (a shorter time than the results without the ramp limit). 


VII. Conclusions 

This paper has presented a discussion of the tuning techniques used for the current regulator of the NASA GRC 
flywheel system. The flywheel system electronics incorporate an AC filter between the inverter and the 
motor/generator to reduce the dv/dt stress on the motor and to improve the magnetic bearing control. The AC filter 
adds additional dynamics to the system that must be considered when tuning the current regulator in order to achieve 
the desired performance at high speeds. The current regulator was initially tuned without consideration of the AC 
filter and the results were satisfactory at low speeds. However, as the speed increased, the response to a step change 
in current command was found to be more and more oscillatory. The oscillations were eliminated and a fast 
response was achieved when the current regulator gains were changed to include the impedance of the AC filter. The 
response with the adjusted gains was improved further by augmenting the current regulator with a back EMF 
decoupling portion and limiting the rate of change of the commanded current. A step change in current command is 
not a realistic command trajectory and limiting the current command to a ramp achieves the desired bandwidth with 
less overshoot and settling time. The ramped command also reduces the voltage command overshoot during the 
transient thus ensuring that the commanded voltage can be achieved at all times. 

A few discrepancies were noted that were most likely caused by an error in the estimated rotor position. Adding 
back EMF decoupling to the current regulator did not reduce the transient overshoot in the actual phase current even 
though it was reduced in the q- and d-axes. For a perfect transformation, the phase current transient would have also 
been reduced. Additionally, the commanded voltage was lower than the back EMF voltage at high speeds during 
acceleration. Again, with perfect transformation using the exact rotor position angle, the commanded voltage is 
expected to be higher than the back EMF voltage during acceleration. Quantifying and correcting the suspected rotor 
position estimation error is the next step for the improvement of the flywheel motor/generator control. 

However, the overall goal of achieving current regulation with a 2 kHz bandwidth for speeds up to 50,000 RPM 
was achieved and documented. This type of current response enables the high performance torque control necessary 
for the flywheel operation at all speeds. 
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Figure 31. Commanded Phase Voltage During Step 
Change in Current Command with 60 kA/sec Rate 
Limit. 



time, seconds 

Figure 32. Motor Phase Current with 60 kA/sec Rate 
Limit at 50,000 RPM. 
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Figure 33. Motor q-axis Current with 60 kA/sec Rate 
Limit. 



Figure 35. Motor Speed During Step Change in 
Current Command with 60 kA/sec Rate Limit. 



time, seconds 

Figure 34. Magnified View of Motor q-axis Current 
with 60 kA/sec Rate Limit. 



Figure 36. Magnified View of Motor d-axis Current 
with 60 kA/sec Rate Limit.. 
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Appendix 

The Thevenin equivalent circuit was calculated for the filter shown in Fig. 3. The results are given below where 
Vthev, numerator is the Thevenin voltage numerator, V t hev, den o m i nat o r is the Thevenin voltage denominator, Z, hev>numerator is the 
Thevenin impedance numerator, and Z, h ev, denominator is the Thevenin impedance demoninator. 

Vthev, numerator = S 3 R 2 C 2 L,C t + S 2 L,C t + SR 2 C 2 + 1 (23) 


V thev, denominator 


= s 6 CiC 2 L 1 L 2 C,L, + s 5 LiC 1 R 2 C 2 L,C, + s^CiQL^ + L 2 C 2 L,C, + L,C 2 L t C, + LiCiL,C,) + 
s 3 (L!CiR 2 C 2 + R 2 C 2 L,C,) + s 2 (L,C, + L 2 C 2 + LtC 2 + LA + L,C 2 ) + s(R 2 C 2 ) F 1 


(24) 


z, hev , numerator = s (C t C 2 C,L t L 2 L,R 2 ) + s (C 1 C,L 1 L 2 L,) + S 4 (C 2 C,L 1 L,R 2 + C 2 C,L 2 L,R 2 + C!C 2 L 2 L,R 2 + C 1 C 2 L 1 L 2 R 2 ) + 
s 3 (C,L 1 L, + C,L 2 L, + CiL 2 L, + CiL[L 2 ) f s 2 (C 2 L,R 2 +C 2 LiR 2 + C 2 L 2 R 2 ) + s(Lt+L 2 +L,) (25) 


Zthev, denominator S (CtC 2 CtL[L 2 Lt) F S (L | L?C [L] L t R 2 ) F S (C 2 C)Li L[ F C | C 2 L ] L 2 + C | C [L| L t F C 2 CtL 2 Lt F CtC 2 L 2 Lt) F 

s 3 (CiC 2 L,R 2 F C 2 C,L,R 2 F CiC 2 LiR 2 ) F s 2 (L,C 2 F L,C, F CiL, F C , L, F C 2 L 2 FC 2 Li) F sC 2 R 2 F 1 (26) 


For the filter without the L-C trap, the Thevenin values are as follows. 


= sR 2 C 2 f 1 

V thev.no L-C trap + S 3 (L 1 C 1 R 2 C 2 ) F S 2 (L 2 C 2 F LtC 2 F LA) F S (R 2 C 2 ) F 1 

s 4 (C 2 CiLiL 2 R 2 ) f s 3 (CiLiL 2 ) f s 2 (C 2 LiR 2 f C 2 L 2 R 2 ) f s(LiFL 2 ) 
thev.no L-C trap " + S 3 (C ALiR 2 ) F S^CA F C 2 L 2 FC 2 Li) F S C 2 R 2 F 1 


(27) 

(28) 
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